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Abstract 

Vortices in supersymmetric gauge field theory are important constructs in 
a basic conceptual phenomenon commonly referred to as the dual Meissner 
effect which is responsible for color confinement. Based on a direct minimiza- 
tion approach, we present a series of sharp existence and uniqueness theorems 
for the solutions of some non-Abelian vortex equations governing color-charged 
multiply distributed flux tubes, which provide an essential mechanism for lin- 
ear confinement. Over a doubly periodic domain, existence results are obtained 
under explicitly stated necessary and sufficient conditions that relate the size of 
the domain, the vortex numbers, and the underlying physical coupling parame- 
ters of the models. Over the full plane, existence results are valid for arbitrary 
vortex numbers and coupling parameters. In all cases, solutions are unique. 



1 Introduction 

A fundamental puzzle in physics, known as the quark confinement, is that quarks, 
which make up elementary particles such as mesons and baryons, cannot be observed 
in isolation. A well accepted confinement mechanism, known as the linear confinement 
model, states that, when one tries to separate a pair of quarks, such as a quark and 
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an anti-quark constituting a meson, the energy consumed would grow linearly with 
respect to the the separation distance between the quarks so that it would require an 
infinite amount of energy in order to split the pair. The quark and anti-quark may be 
regarded as a pair of source and sink of color-charged force fields. The source and sink 
interact through color-charged fluxes which are screened in the bulk of space but form 
thin tubes in the form of color-charged vortex-lines so that the strength of the force 
remains constant over arbitrary distance, resulting in a linear dependence relation 
for the potential energy with regard to the separation distance. Such a situation 
is similar to that of a magnetic monopole and anti-monopole pair immersed in a 
type-II superconductor. The magnetic fluxes mediating the interacting monopoles 
are not governed by the Maxwell equations, which would otherwise give rise to an 
inverse-square-power law type of decay of the forces and lead to non-confinement, but 
rather by the Ginzburg-Landau equations, which produce thin vortex-lines, known as 
the Abrikosov vortices or the Nielsen-Olesen strings. The repulsion of the magnetic 
field in the bulk region of the superconductor is due to the Meissner effect and the 
partial magnetic penetration of the superconductor in the form of vortex-lines is a 
consequence of type-II superconductivity. Thus, one may conclude that a monopole 
and anti-monopole immersed in a type-II superconductor would be linearly confined. 

Inspired by the above- described monopole confinement in a type-II superconduc- 
tor, Mandelstam [36, 37], Nambu [39], and 't Hooft [54, 56] proposed in the 1970s that 
the ground state of quantum-chromodynamics (QCD) is a condensate of chromomag- 
netic (color-charged) monopoles, causing the chromoelectric fluxes between quarks to 
be squeezed into narrowly formed tubes or vortex-lines, similar to the electron con- 
densation in the bulk of a superconductor, in the form of the Cooper pairs, resulting 
in the formation of magnetic flux-tubes or vortex-lines which mediate the interaction 
between monopoles, following a non-Abelian version of the Meissner effect, called 
the 'dual Meissner effect' [9, 20, 33, 52], which is responsible for the screening of 
chromoelectric fluxes [49, 50]. 

Interestingly, although finite-mass monopoles of the 't Hooft [53] and Polyakov 
[43] type in non-Abelian gauge theory have long been demonstrated to exist [2, 18], 
the magnetic fluxes are Coulomb-like which spread out radially and will not give rise 
to confinement. In 1994, Seiberg and Witten [45] came up with an M = 2 supersym- 
metric gauge-field-theoretical formalism of non-Abelian monopole condensation and 
studied its implications to color confinement. Later, in 1997, Hanany, Strassler, and 
Zaffaroni [25] showed that the flux tubes or strings produced in the Seiberg- Witten 
formalism are of the Abrikosov-Nielsen-Olesen (Abelian) type [1, 23, 40, 57, 58] 
which are not exactly what anticipated in QCD [25, 50, 49]. In 2002, Marshakov 
and Yung [38] constructed Abelian vortices in a softly broken M = 2 supersymmet- 
ric QCD (SQCD) model and showed that, although confinement is due to Abelian 
flux tubes, the multiplicity of the meson spectrum is the same as expected in a the- 
ory with non-Abelian confinement. In 2003, Hanany and Tong [26] derived a broad 
class of non-Abelian vortex equations and computed the dimensions of the associ- 
ated moduli spaces, and Auzzi, Bolognesi, Evslin, Konishi, and Yung [5] analyzed 
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non-Abelian vortices and confinement in M = 2 SQCD in the context of non-Abelian 
superconductors. Since then, the subject of non-Abelian vortices, monopole conden- 
sation, and confinement has been extensively developed [12, 13, 15, 22, 27, 47, 48]. 
See [14, 21, 32, 46, 49, 50] for surveys and further literature. Mathematically, these 
studies unveil a broad spectrum of systems of elliptic equations with exponential 
nonlinearities and rich properties and structures, which present new challenges. 

Recently, C. S. Lin and one of the authors (Y. Y.) carried out a systematic study 
[34, 35] of the multiple vortex equations obtained in [13, 15, 14, 22, 47, 48, 49, 50]. 
A series of sharp existence and uniqueness theorems were established. The methods 
used include monotone iterations, a priori estimates and degree-theory argument, and 
constrained minimization. In the present paper, we do two things. One is to develop 
and prove another series of sharp existence and uniqueness theorems for the multiple 
vortex equations derived in [5, 6], which are not covered in [34, 35]. The second thing is 
to develop a methodology that has not normally been used for these kinds of problems, 
over doubly periodic domains, which are often more difficult to approach due to the 
appearance of some integral constraints naturally associated with the equations. This 
is the highly efficient direct minimization approach which enables us to identify the 
key analytic ingredients and pursue a complete understanding of the problems almost 
immediately. As a by-product, such an approach also provides a constructive method 
for solutions. It is hopeful that our method may be explored further to study various 
multiple vortex equations, arising in non-Abelian gauge field theory, of more difficult 
structures. It should be noted that, after solving the problems here by the direct 
method, we gained true insight to solve them by the usual constrained minimization 
method. But we were not successful in all cases. Thus, the direct method, besides 
being simpler, is sometimes the only constructively workable one. 

The content of the rest of the paper is outlined as follows. 

In Section 2, we recall the SQCD multiple vortex equations of Auzzi, Bolognesi, 
Evslin, Konishi, and Yung [5] obtained in 2003 in which vortices are induced from 
three complex scalar fields. We then state our existence and uniqueness theorem, In 
the next two sections, we give the proofs for the various parts of the theorem. In 
Section 5, we turn our attention to a study of the multiple vortex equations derived 
by Auzzi and Kumar [6] in a supersymmetric Chern-Simons-Higgs theory formulated 
by Aharony, Bergman, Jafferis and Maldacena [3] , known as the AB JM model. In this 
problem, vortices are generated from m complex scalar fields and the governing elliptic 
system consists of m equations. Again, we are able to obtain a sharp existence and 
uniqueness theorem for solutions over a doubly periodic domain and the full plane. 
Proofs of results are sketched in the subsequent two sections, as earlier. 

In Section 8, we further illustrate how our direct methods may be used in tackling 
other problems of similar structures. Specifically, in §8.1, we revisit the SO(2N) BPS 
vortex equations of Gudnason, Jiang, and Konishi [22] , studied in [35] where solutions 
are constructed by a constrained minimization method when the total vortex number 
n does not exceed 3 and existence of solutions is established for arbitrary n using a 
degree-theory argument. 
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The main difficulty encountered in [35] is an extra term in one of the equations 
that makes it hard to resolve the constraints explicitly, which may be flown away by 
a homotopy flow. Here is an example where the direct method seems to be the only 
(constructively) successful one. With it we are able to obtain an existence proof for 
an arbitrary n. In fact, we will carry out our study in the most general situation 
where vortices are induced from the sets of zeros of the two complex scalar fields of 
the model. 

In §8.2, we present a sharp existence and uniqueness theorem for the multiple 
vortex equations obtained by Marshakov and Yung [38] in 2002, which may be re- 
garded as the earliest non-Abelian SQCD vortex equations for which the vortex-lines 
are taking values in the Cartan subalgebra of £77(3) and, also, the starting point of 
the later development of the subject of non-Abelian vortices and monopoles in SQCD 
and their applications to color confinement. The method is again centered around 
direct minimization. 

2 Vortices in Yang— Mills— Higgs theory 

Following Auzzi, Bolognesi, Evslin, Konishi, and Yung [5], the Yang-Mills-Higgs 
action hosting the gauge field theory undergoing the spontaneous symmetry breaking 

SU(N) ->• SU(N- 1) x 1/(1), (2.1) 

within the context of the critical BPS coupling, assumes the form 

s = /{-^^-^ + (v,?a) , vV 

-y G&V) 2 - m K + 1) (^ A - m 2 } dx, (2.2) 

where K — N — 1, the index a = 1, 2, • • • , K 2 — 1, labels the group generators {t a } of 
SU(K), g, e > are the SU(K) and U(l) gauge-field coupling constants, respectively, 
£ > determines the vacuum expectation value of the quark fields q A lying in the 
fundamental representation of SU(K) x U(\), A — 1, 2, • • • , iV flavor runs over the quark 
flavors, 

v„ = a,. - up - iV, <° - 7== ( \ K - K ) • < 2 - 3 > 

denotes the gauge- covariant derivative, with l m the m x m identity matrix, and 

the Minkowski spacetime is of the signature (H ). As a consequence of the 

Bogomol'nyi reduction [8, 31] for static vortex solutions, the BPS vortex equations 
[5] are of the form 

\ F ? 2 + (q\t a q A ) = 0, a = l,2,-..,X 2 -l, (2.4) 
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- 2 F l2 + {q\q A -Ki) = 0, (2.5) 



e 



V iq A + iV 2 q A = 0, A = l,2,---,iV flavor . (2.6) 



In the specific situation of N = 4 so that the unbroken symmetry is given by the 
group SU(3) x i7(l), the non-Abelian vortex solutions may be described by gauge 
fields solely given in the a = 3, 8 (these are the Cartan subalgebra indices in the 
Gell-Mann matrix representation) and the {7(1) sectors, and the quark fields are 
represented by the complex matrix 

(q kA )= I r I . (2.7) 




where k = 1,2,3, is the color index which runs vertically, A = 1,2,3, is the flavor 
index which runs horizontally, and the winding numbers of 0, ip, x, away from a local 
region where (j>,tp,x ma Y vanish, say rii,n 2 ,n 3 , characterize the quark fields, which 
will be identified as the vortex charges or vortex numbers. 
Now set 

A) = a v A* = ^=b 3 , A 3 = l -c 3 , j = 1,2. (2.8) 

where dj,bj,Cj (j = 1,2) are real- valued vector fields. Then, in terms of the fields 
given in (2.7) and (2.8), the non-Abelian BPS multiple vortex equations (2.4)-(2.6) 
are found [5] to be 

(d 1 + id 2 )(f) = i Q(ai + ia 2 ) + i(6i + ib 2 ) + i( Cl + ic 2 )^) 0, (2.9) 

(9x + i9 2 )^ = i^-i(a 1 + ia 2 ) + i(6i + i6 2 ) + ^(ci + ic 2 )^, (2.10) 

(d 1 + id 2 ) X = i(^-^(6 1 + i6 2 ) + i(c 1 + ic 2 )^x, (2.H) 

rv 

ai2 = "2(I0| 2 -H 2 ), (2-12) 

h2 = --(|0| 2 + H 2 -2| X | 2 ), (2.13) 

ci2 = -/3(|0| 2 + I^| 2 + |X| 2 -3O, (2.14) 

where 

ai2 = d\a 2 — d 2 a±, etc, (2-15) 

are the reduced field curvatures, and a = g 2 ,(3 = 3e 2 . For convenience, we now use 
the complexified variables 

z = x 1 + ix 2 , a = ai + ia 2 , b = b x + ib 2 , c = C\ + ic 2 , (2-16) 
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and the complex derivatives 



d = d z = ±(d 1 -id 2 ), d = d g =±(d 1 +id 2 ) (2.17) 
to convert (2.9)-(2.14) into the system 

^ = I(-| + J + £), (,19) 

3lnX = i(-M), P.20) 

i(da-9a) = |(|0| 2 -|Vf), (2.21) 

i(db-M) = |«| 2 +|</f-2| X | 2 ), (2.22) 

i(Sc-9e) = /3(H 2 + h/>| 2 +|x| 2 -3£), (2.23) 

away from the zeros of (f>,ip,x- From (2.18)-(2.20), we have 

a = -20(ln0-lnV), (2.24) 

b = -20(ln0 + ln^-21nx), (2.25) 

c = -20(ln0 + ln^ + lnx). (2.26) 

Inserting (2.24)-(2.26) into (2.21)-(2.23) and using the result A = Add = Add, we 
obtain 

A(ln|0| 2 -lnH 2 ) = a (|0| 2 -|^| 2 ), (2.27) 

A(ln|0| 2 + ln|^| 2 -21n| X | 2 ) = a(|0| 2 + |^| 2 - 2| X | 2 ), (2.28) 

A(ln|0| 2 + ln|^| 2 + ln| X | 2 ) = 2/3(|0| 2 + |Vf + \x\ 2 - 30, (2-29) 

again away from the zeros of <p,if},x- Following [31], we know that the equations 

(2.9)— (2.1 1) or (2.18)-(2.20) imply that the zeros of 0, ip, x are discrete and of integer 
multiplicities. We use Z^, Z^, Z x to denote the sets of zeros of 0, ip, x, 

Z4, = {pi,i,---,Pi,m}, Z^ = {p 2 ,i,---,p 2 ,n 2 }, Z x = {p 3) u- • • >P3,n 3 }, (2-30) 

so that the repetitions among the points pe tS , £ — 1, 2, 3, s — 1, • • • , ri£, take account 
of the multiplicities of these zeros. Then the substitutions 

Mi = ln|0| 2 , n 2 = In |V| 2 , u 3 = In \ X \ 2 , (2.31) 

enable us to recast the equations (2.27)-(2.29) into the following elliptic system 

X>i>)-X>»>) , (2-32) 

5 = 1 S=l / 



A(ui + u 2 - 2u 3 ) = a(e ui + e" 2 - 2e U3 ) 

E <w (*) + E <w w - 2 E <w (*) '( 2 - 33 ) 

s=l s=l s=l / 

A( Mi + M2 + m 3 ) = 2/3(e Ul + e na + e" a - 30 

E <W (*) + E <W (*) + E <W (*) , (2.34) 

s=l s=l s=l / 

now denned over the entire domain. 

Two situations are of interest, namely, the situation where the equations are con- 
sidered over a doubly periodic domain, Q, governing multiple vortices hosted in fl 
so that the field configurations are subject to the 't Hooft periodic boundary condi- 
tion [55, 60, 61] under which periodicity is achieved modulo gauge transformations, 
and the situation where the equations are considered over the full plane M 2 and the 
solutions satisfy the boundary condition 

u e (x)->]n£ as \x\ ->■ oo, £=1,2,3. (2.35) 

Concerning these situations, our main existence and uniqueness theorem for solu- 
tions of (2.9)-(2.14) or (2.32)-(2.34) may be stated as follows. 

Theorem 2.1. Consider the BPS system of multiple vortex equations (2.9)-(2.14) 
for (<f),ifj,x,cij,bj,Cj) with the prescribed sets of zeros given in (2.30) so that <p,ifj,x 
have ni,n 2 ,n 3 arbitrarily distributed zeros, respectively. 

(i) For this problem over a doubly periodic domain fl, a solution exists if and only 
if the following three conditions 

^-(ni + n 2 - 2n 3 ) + -(m - n 2 ) + -3(^1 + n 2 + n 3 ) < (2.36) 
6a a 6p 2n 

1 . . 1. i . . 

— (ni+n 2 -2n 3 ) (ni - n 2 ) + -^(ni + n 2 + n 3 ) < ——, (2.37) 

6a a 6p 2ir 

_A( ni+n2 _2n 3 ) + ^r(n 1 +n 2 + n 3 ) < ^p, (2.38) 
6a 6p 2ir 

hold simultaneously. Moreover, whenever a solution exists, it is unique. 

(ii) For this problem over the full plane M. 2 subject to the boundary condition 

I0| 2 ,M 2 ,|X| 2 ^£ as \x\^oo, (2.39) 

there exists a unique solution up to gauge transformations so that the boundary be- 
havior stated above is realized exponentially rapidly. 

In either case, the excited total vortex fluxes are quantized quantities given explic- 
itly by the formulas 

$«= [a 12 dx = 27r(m-n 2 ), (2.40) 
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= Jb 12 dx = 2vr(ni +n 2 - 2n 3 ), (2.41) 
$ c = c 12 dx = 27r(ni + n 2 + n 3 ), (2.42) 



respectively 

This theorem will be established in the following two sections. 

3 Proof of existence for doubly periodic case 

In this section, we consider the equations (2.32)-(2.34) defined over a doubly periodic 
domain, Q. Let w° be a solution of the equation 

A^ = -^ + 47r^^>), xeQ, £ = 1,2,3. (3.1) 



s=l 



Then the substitutions 

■u^ = w° + i^, £ = 1,2, 3, 
recast the equations (2.32)-(2.34) into 

A( Vl -v 2 ) = a(e u ° +n -e u ° + ^) + p(m-n 2 ), (3.2) 

A( Vl + t; 2 - 2v 3 ) = a(e u ° +Vl +e u ° +V2 -2e u ° +V3 ) + ^(n 1 + n 2 -2n 3 ), (3.3) 

A(ui +v 2 + v 3 ) = 2P(e< +vl +e u ° +v *+e u ° +V3 -30 + ^(n 1 +n 2 + n 3 ).(3A) 

Naturally, we should use the transformation 

= v 1 - V 2 , ( Vi = \wi + \w 2 + \w z , 

= vi+v 2 - 2v 3 , < v 2 = -\wi + hv 2 + §w 3 , (3.5) 
= Vi + v 2 + v 3 , [ v 3 = -|iu 2 + 3W3, 

to change (3.2)-(3.4) into the equations 

Aw! = a {e u ° + ^ wl+ ^ W2+ ^ m - e u °-^ wl+ ^ W2+ ^ ws ) + ^(m -n 2 ), (3.6) 
At« 2 = o;(e"i >+ ^ 1+ l W2+ l W3 + e M 2"5 wl+ l W2+ l W3 - 2e u 3-| W2+ i"' 3 ) 

47T 

+p(ni + n 2 -2n 3 ), (3.7) 

Au> 3 = 2j3(e u ° + ^ Wl+ ^ W2+ ^ m + e u 2"i wl+ ^ 2+ 5"' 3 + e u *~* W2+ * W3 - 3£) 

+|^|(ni + n 2 + n 3 ), (3.8) 




which are easily seen to be the Euler-Lagrange equations of the functional 

I(w 1: w 2 ,w 3 ) = [ (^|V Wl | 2 + ^|V W2 | 2 + ^-|V W3 r + e M " + ^ 1+ ^ 2+ ^ 3 
J n {4a 12a 12(3 

+( A-± wl +\w 2 +\w 3 + e «§-i»2+i»3 + n r ni _ n2 ) Wl 

a\Q\ 

+ ^k {ni +n2 - 2n3)W2 + {wk {ni +n2+ ns) - e ) W3 } dx - 

(3.9) 

This functional is not bounded from below when 

27r(m + ra 2 + n 3 ) > 3/3£|ft|. (3.10) 

In fact, we can show that (3.10) will never happen for (3.6)-(3.8). Indeed, integrating 
(3.6)-(3.8), we obtain the conditions 

/ e «i + 3Wl + sW2 + 3W 3 ^ x 

Jn 

( 1 1 1 

= f |f2| - 2tt — (m + n 2 - 2n 3 ) + -(m - n 2 ) + + n 2 + n 3 ) 

\3o; a 3p 

= 77i>0, (3.11) 

( 1 1 1 

= f ft - 2tt — (m + n 2 - 2n 3 ) (m - n 2 ) + — (^i + n 2 + n 3 ) 

\6a a 6p 

V2 > 0, (3.12) 

eU o_i W2+ i W3 ^ 



/ 



/ 2 1 

= £|ft| - 2tt ( + n 2 - 2n 3 ) + g^("i + n 2 + n 3 ) 

= %>0, (3.13) 

which are exactly the conditions (2.36)-(2.38). In particular, we have 

^1+^2 + % >0, (3.14) 

which rules out (3.10) immediately. 

Below, we shall show that, under the conditions (3.11)— (3.13), the equations (3.6)- 
(3.8) have a solution. We will use both a direct minimization method and a con- 
strained minimization method to approach the problem. These methods may be of 
independent practical value for computational purposes. 
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3.1 Direct minimization 

We use W 1,2 (Q) to denote the usual Sobolev space of scalar- valued or vector- valued 
fi-periodic L 2 -functions whose derivatives are also in L 2 (Q). In the scalar case, we 
may decompose W 1>2 (Q) into W 1,2 (ft) = R © W 1,2 (Q) so that any / G W^ 2 (VL) can 
be expressed as 

(3.15) 



/ = / + /, feww(n), / fdx = o. 

Jn 

It is useful to recall the Moser-Trudinger inequality [4, 16] 

( — [ \Vu\ 2 dx] , u e W 1 ' 2 ^). 
V 167r Jn J 



/ e"dx < Cexp 
Jn 



(3.16) 



With (3.16), it is clear that the functional / defined by (3.9) is a (^-functional 
with respect to its argument (wi, w 2 , w 3 ) G W 1,2 (Q) which is strictly convex and lower 
semi-continuous in terms of the weak topology of W l,2 (Q). 

With the notation (3.15), we may apply the transformation (3.5) to arrive at 



I(w 1 ,w 2 ,w 3 ) 



4a 



|Vwi| 



12a 



\Vw 2 \ 2 + 



12/3 



dx 



^1 + ln 


'at 







(3.17) 



where we have used the Jensen inequality to obtain the lower bounds 

J e u° e +vt+vt dx > |Q| exp J^JL J ( u + ^ + ^) dx' 



Q\ exp 



1 

IfTl 



■u* dx 



eiu = ae e^, £= 1,2,3, (3.18) 



in (3.17). Thus, in particular, we see that / is bounded from below and we may 
consider the following direct minimization problem 

r] = inf {l( Wl ,w 2 ,w 3 ) w u w 2 ,w 3 G W ll2 (fi) j . (3.19) 
Let {(w^, w^)} be a minimizing sequence of (3.19). Since the function 



F(u) = ae u - i]u, 



(3.20) 



where a, i] > are constants, enjoys the property that F(u) — > 00 as u — > ±00, we 
see from (3.17) that the sequences (£ = 1,2,3), hence {y^} {£ = 1,2,3), are 

all bounded. Without loss of generality, we may assume 



wi"- 1 — > some point u/f°^ G R as n — > 00, £ = 1, 2, 3. 



(3.21) 
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On the other hand, in view of (3.17) and the Poincare inequality, we see that all 
the sequences {wf 1 ^} are bounded in W 1,2 (Q), £ = 1,2,3. Without loss of generality, 
we may assume 

->■ some element G W h2 {tt) weakly as n ->■ oo, £=1,2, 3. (3.22) 

Of course, G W 1 ' 2 ^) {£=1,2,3). Set =u/f o) +w < f o) {£ = 1,2,3). Then 
(3.21) and (3.22) lead us to ->• ti/ 00 ) weakly in W^ft) as n ->■ oo {£ = 1,2,3). 
The weakly lower semi-continuity of / enables us to conclude that {w[°°\ w^, w^) 
solves (3.19), which is a critical point of /. As a critical point of I, it satisfies the 
equations (3.6)-(3.8). Since I is strictly convex, it can have at most one critical point. 
Thus, the uniqueness of the solution of (3.6)-(3.8) follows immediately. 



3.2 Constrained minimization 

For convenience, we rewrite the constraints (3.11)— (3.13) collectively as 

J E {w 1 ,w 2 ,w 3 ) = [ e u ° +v * dx = r] e , £ = 1,2,3, (3.23) 
Jn 

and consider the constrained minimization problem 

r]o ee inf ^I{w 1 ,w 2 ,w 3 ) {w 1 ,w 2 ,w 3 ) G W 1 ' 2 ^) and satisfies (3.23) | . (3.24) 

Suppose that (3.24) allows a solution, say (wi, w 2 , w 3 ). Then there are numbers 
(the Lagrange multipliers) in R, say Ai, A2, A3, such that 

D(I + A1J1 + \ 2 J 2 + \ 3 J 3 ){w 1 ,w 2 ,w 3 ){w 1 ,w 2 ,w 3 ) = 0, 

y{w\,w 2 ,w 3 ) G W 1)2 {VL). (3.25) 

Now take the trial configurations, {w\,w 2 ,w 3 ) = (1, 0, 0), (0, 1, 0), (0, 0, 1) consecu- 
tively, in (3.25). Since {w\, w 2 , w 3 ) satisfies the constraints (3.23), as a result, we 
have 

1 1 1 1, 1 1 1, 1, 

2 A i^i _ 2 2X12 = ' 6 1 + 6 2 ~~ 3 3 = 3 1 + 3 2 + 3 3 = 

(3.26) 

Consequently, Ai = A 2 = A 3 = 0. In other words, the Lagrange multipliers disappear 
automatically and, thus, the search for a critical point of the functional I is converted 
to obtaining a solution of the constrained minimization problem (3.24). 
In order to approach (3.24), we resolve (3.23) to write down 



2k 



hir^-ln (^J e u ° +ie dx^J , £=1,2,3. (3.27) 
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Hence, in view of the left-hand side of (3.17), we get 

> - Ve ^ m + ^2 m In ( f e u ^ dx\ > ^ m In f^j , (3.28) 
e=i i=i ' e=i We/ 

where we have used the Jensen inequality and the definition of the quantities <Te 
(£ = 1,2,3) given in (3.18). Thus the problem (3.24) is well defined. 

Let {(w^Wj™',^™')} be a minimizing sequence of (3.24). Then (3.28) says that 
the sequence {(w^\w^ l \w^)} is bounded in W 1,2 {Vt). Hence we may assume that 
{(w^, w^jtb^)} is weakly convergent in W 1,2 (Q). The inequality (3.16) and the ex- 
pressions (3.27) indicate that {(w^\ w^ 1 ^)} is a convergent sequence in M 3 . Thus 
{(w^, )} is weakly convergent in W 1,2 (Q). In view of the weak continuity 

of the constraint functionals Jg defined in (3.23) and the weak lower semi-continuity 
of the functional I defined in (4.11), we see that the weak limit of {(w^, w^, Wg^)} 
in W l,2 (Q) is a solution of (3.24). As a critical of J, it is also unique. Therefore, a 
constrained minimization proof for the existence of a unique solution of the equations 
(3.6)-(3.8) is obtained. 



4 Proof of existence for planar case 

With the correspondence relations stated in (3.5), we have 

3wl + wl + 2wj = Q{vl + vl + vl), (4.1) 
3|V Wl | 2 + |V«; 2 | 2 + 2|V«;3| 2 = 6(|V^| 2 + |V^ 2 | 2 + |V^ 3 | 2 ), (4.2) 

which will be useful for our analysis to follow. 

To proceed further, here and elsewhere in the paper when we deal with the planar 
cases, we use the method developed in [31] and introduce the background functions 
[31] 

s=l 

(Here and in the sequel, the parameter /i > should not be confused with the 
spacetime index \i used in the context of various field equations.) 
Then we have 

ri£ rig 

Au? = -/* + 47r5%», h e (x)=4j2 £ = 1,2,3. (4.4) 

s=1 s=1 \P -r \ x Pt,s\ ) 

Now use the substitutions 

u e = \n£ + u° e +v e , £=1,2,3, (4.5) 
a£ ^ a, P^P, (4.6) 



= -^ln(l + fi\x -p £tS \~ 2 ), fi>0, £ = 1,2,3. (4.3) 
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Aw 2 



and (3.5) in (2.32)-(2.34), we obtain the governing equations 

Aw! = a (e^i+h^+k^H^ - e u° 2 -^ 1+ i W2+ i W3 ^ + ^ _ ^ ^ 

- q, ^ e «l + 5Wl + g^2 + |«)3 _|_ e U 2 -\ Wl + lw2 + \w 3 _ 2 e «3-3 w 2 + 3W 3 ^ 

+(/ii + /i 2 -2/i 3 ), (4.8) 

+(h 1 + h 2 + h 3 ), (4.9) 

over the full plane R 2 . The boundary condition for Wi,w 2 ,w 3 reads 

w e (x)->Q as|x|->oo, £=1,2,3. (4.10) 

In order to obtain a solution of (4.7)-(4.9) subject to the boundary condition 
(4.10), we look for a critical point of the action functional 



-Wi H — w 2 H — w 3 
2 6 3 

1 1 1 

+ -w 2 + -w 3 

+ 77-^1 - ^2)^1 + 7— (^1 + ^2 - 2/i 3 )w 2 + T^(/ii + h 2 + h 3 )w 3 \ dx, 
la ba bp J 

(4.11) 

which is C 1 and strictly convex over W^ 1,2 (R 2 ). After some algebra, it can be seen 
that the Frechet derivative of / enjoys the following property, 



(DI(w 1 ,w 2 ,w 3 ))(w 1 ,w 2 ,w 3 ) - 



— \Vw 1 \ 2 + — \Vw 2 \ 2 + 
2a ba 



6/3 



. ( 3 3 3 1 

" /r2 [ t=\ e=i t=\ J 



dx 



(4.12) 



where g/s are some linear combinations of he's. Now setting 7 = maxja, 2/3} and 
applying (4.2), we derive from (4.12) the lower bound 

(DI(w 1 , w 2 , w 3 ))(w u w 2 , w 3 ) 



> 



[ \ - E 1 v ^i 2 + E i eU ^ vl - % + [ e " ? - % + 1 d;r 

/ { - E I V^| 2 + E - 1 + ft) } dx. 



(4.13) 
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We can now follow the analysis in [31]. To simplify the notation, we suppress the 
subscript t and rewrite a typical part on the right-hand side of (4.13) as 



M(v) 



[ v (e ui>+v -1 + g) dx. 
Jr 2 ^ ' 



(4.14) 



Here g should not be confused with the coupling constant used before in the field- 
theoretical context. Thus e* — 1 > t (t E R) gives us 



u u +v 



l + g>u° + v + g, 



which leads to 

M(v + ) > vldx+ v + (u° + g)dx 
Jr 2 Jr 2 

> \ \ v\dx-\\ (u° + g) 2 dx. (4.15) 
z Jr 2 1 Jr 2 

On the other hand, in view of the inequality 1 — e - ' > t/(l + 1) (t > 0), we have 
v„(l- g-e u °~ v -^j = v. (l-g + e u °[l-e- v -] - e M °) 



> v-[l-g + & 



,0 V- 



1 + v. 



TTir( 1 -») + Tfir( 1 - rf '- 9 )- (4 ' 16) 



Of course, we may choose /i > in (4.3) large enough so that g < 1/2 (say). Further- 
more, since 1 — e u and g are in L 2 (R 2 ), we have 









l-e u °-g 


dx<-l 


/ R 2 1 + V_ 


4 h 



Combining (4.16) and (4.17), we obtain 

M( -"- ) = I v ~ i 1 ~ 9 " e "°~*~) dx £ 1 1 TTvI dx ~ c ' (418) 

where and in the sequel C > denotes an irrelevant constant. Summarizing (4.15) 
and (4.18), we arrive at 



M(v) > 



1 



dx-C. 



4 Jr2 1 + \v\ 
We now recall the standard Sobolev inequality 

v A dx<2 v 2 dx \Vv\ 2 dx, v E W 



1,2/Tn>2\ 



(4.19) 



(4.20) 
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Consequently, we have 

/ v 2 dx I = 

Jr 2 J 

< 2 

< 4 



1 

< - 
~ 2 



\v\ 


('I 


1 + \v\ K 
v 2 


(l+ f 

V' 


> 


(1 + 





\v\)\v\ dx 



dx 



[ (v 2 + v 4 ) 
Jr 2 



dx 



dx I v dx I 1 + 

>2 



da; + C 1 + 



[i + \v\y 



\Vv\ 2 dx 



dx 



+ 



|Vi>| dx 



(4.21) 



\Vv\ 2 dx + 



Jr 2 (1 



dx 



(4.22) 



As a result of (4.21), we have 

(f v 2 dx) 2 <c(l+ _ , 

Now set C = min{l/7, 1/4}. In view of (4.13) and (4.19), we have 

(DI(w 1 ,w 2 ,w 3 ))(w 1 ,w 2 ,w 3 )>C y] [ (W| 2 + -^— ) dx-C. (4.23) 

JR 2 \ 1 + 1^1 / 

As a consequence of (4.22), (4.23), (4.1), (4.2), we may conclude with the coercive 
lower bound 



(DI(w 1 ,w 2 ,w 3 ))(w 1 ,w 2 ,w 3 ) > CiJ^ H^Hh/1,2 



C 2 , 



(4.24) 



where C 1: C 2 > are constants. In view of the estimate (4.24), the existence of a 
critical point of the functional / in the space W /1 ' 2 (IR 2 ) follows. In fact, from (4.24), 
we may choose R > large enough so that 

inf i^(DI(w 1 ,w 2 ,w 3 ))(w 1 ,w 2 ,w 3 ) | ^ H^Ivf 1 - 2 ^ 2 ) = -R j > 1 (4.25) 

(say). Consider the minimization problem 

3 

(4.26) 



i)o = inf <| I(w 1 ,w 2 ,w 3 ) | ^ ||^|| w i,2(R2) < R 



This problem obviously has a solution due to the fact that the functional (4.11) is 
weakly lower semi-continuous. Let (wi,w 2 ,w 3 ) be a solution of (4.26). We show that 
it must be an interior point. Otherwise, if 



e=i 



R, 



(4.27) 
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then, with the vector notation w = (wi, w 2 , w 3 ), the result (4.25) gives us 

^ J([l-*]w)-J(w) = _ = _ (£)JW)W < _ L (4 . 28) 

Thus, when t > is sufficiently small, with w* = (1 — t)w, we have 

3 

I(w\,wl,wl) < I{w u w 2 ,w 3 ) = rj , ^ 11^11^1,2^2) = (1 - t)R < R, (4.29) 

i=\ 

which contradicts the definition of i] made in (4.26). Thus (wi, w 2 , w 3 ) must an inte- 
rior point for the problem (4.26). Consequently, it is a critical point of the functional 
(4.11). The strict convexity of the functional implies that such a critical point must 
be unique. In the following, we rewrite wg as wg. 
Besides, using the standard embedding inequality 

p-2 

||/|U W < H - l]) 2P \\f\\w^), p > 2, (4.30) 
and the MacLaurin series 



oo 



- 1) 2 = / 2 + E ^z*. ( 4 - 31 ) 

s=3 S - 

it is seen that e / - 1 G L 2 (M 2 ) when / G W 1 ' 2 ^ 2 ). Applying this in (4.7)-(4.9) and 
using elliptic estimates, we have G W^ 2 ' 2 (1R 2 ) (£ = 1, 2, 3). In particular, we(x) — >■ 
as |x| — >■ oo, £ = 1,2,3. In view of this property and (4.7)-(4.9), we see that the 
right-hand sides of (4.7)-(4.9) all lie in L P (R 2 ) for any p > 2, which establishes 
w e G W 2 ' P (R 2 ) (£ = 1,2,3) by elliptic L p -estimates. Consequently, \Vw e \(x) ^ as 
\x\ — > oo {£ — 1,2,3). Linearizing (4.7)-(4.9), we see that wg vanishes exponentially 
fast and Vwe vanishes like O(|o;|~ 3 ) at infinity, I = 1,2,3. Thus, we have 



Aw e dx = 0, £=1,2,3. (4.32) 



Integrating (4.7)-(4.9) over R 2 and inserting (4.32) and the definitions of hg (£ = 
1, 2, 3), we have 

a / ( e Ul + \wi + \w2 + \w3 _ e ul-\wi + \w2 + \w 3 \ fa 

JR 2 ^ ' 

= - [ (h 1 -h 2 )dx = -Air{n 1 -n 2 ), (4.33) 

JR 2 

a / f e U l + h W1 + k W2 + k W3 + qU 2 -\ W1 + ^W2+\w3 _ 2Q u l-\ W2 + \ w, i\ <\ x 
Jr 2 ^ ' 

= - (h 1 + h 2 - 2h 3 ) dx = -47r(ni + n 2 - 2n 3 ), (4.34) 

JR. 2 

2(3 / f e ^i+k w ^+k w ^H w 3 + e^-^i+s^+l^s _|_ gwi-1^2+1^3 _ 3^ ^ 

= -/ (hi + h 2 + h 3 ) dx = -47r(m +n 2 + n 3 ), (4.35) 

JR 2 
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as stated in the theorem. 



5 Vortices in Chern— Simons— Higgs theory 

In the context of supersymmetric Chern-Simons-Higgs theory in the standard (2 + 1)- 
dimensional Minkowski spacetime, recently developed by Aharony, Bergman, Jafferis, 
and Maldacena [3], known also as the ABJM model [6, 7, 19, 24, 51, 59], which is a 
Chern-Simons theory within which the matter fields are four complex scalars, 

C 1 = (Q\Q 2 ,R\R 2 ), I = 1,2,3, 4, (5.1) 

in the bi-fundamental (N, N) representation of the gauge group U(N) x U(N), which 
are all N x N complex matrices, of the gauge fields and B^, and the associated 
Chern-Simons terms for and are set at the levels k and — k so that they give 
rise to the Lagrangian density 

Cos = £^Tr (a^A, + ^A^A, - B,d u B, - ^B,B u B^j , (5.2) 

and the gauge-covariant derivative 

D^C 1 = d^C 1 + iA^C 1 - iC%, 1=1,2,3,4. (5.3) 
The scalar potential density is of the mass-deformed form [19] 

V = Tr(M Qt M a + N a ^N a ), (5.4) 



where 



M a = aQ a + — {2Q^Q\Qft +R?R\Q a -Q a R\R p 



+2Q p R\R a -2R a R\Q p ), (5.5) 

N a = -aR a + —(2R [a R J j 3 R l3] + Q p Q^R a - R a Q^Q p 

+2R?QlQ a -2Q a QlR?), (5.6) 

the Kronecker symbol e a/3 (a, ft = 1,2) is used to lower or raise indices, and a > 

a massive parameter. Thus, when the spacetime metric is of the signature (H ), 

the total (bosonic) Lagrangian density of the ABJM model can be written as 

C = -C cs + TrdD^lD^C 1 ]) - V, (5.7) 

which is of a pure Chern-Simons type for the gauge field sector. As in [6], we focus 
on a reduced situation where (say) R a = 0. Then, by virtue of (5.5) and (5.6), the 
scalar potential density (5.4) takes the form 

V = Tr(M Qt M Q ), M a = aQ a + —{Q a Q\Q p - Q P Q\Q a )- (5.8) 
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The equations of motion of the Lagrangian (5.7) are rather complicated. However, 
in the static limit, Auzzi and Kumar [6] showed that these equations may be reduced 
into the following first-order BPS system of equations 

DqQ 1 - iW 1 = 0, D 1 Q 2 + iD 2 Q 2 = 0, (5.9) 
D X Q X = 0, D 2 Q 1 ^0, D Q 2 = 0, W 2 = 0, (5.10) 

coupled with the Gauss law constraints which are the temporal components of the 
Chern-Simons equations 

A e ^7 F 0B) = i ([ jD Mga]tga -Q^[D»Q a \), (5.12) 

4"7T 



where 



F$ = d,A u -d v A^ + i[A„A u ], 
= d,B„-d u B„ + i[B^B u ], 



W l = aQ 1 + —{Q l Q 2] Q 2 - Q 2 Q 2j Q 1 ), 

K 

W 2 = aQ 2 + ^(Q 2 Q^Qi-Q^Q 2 ), 

K 

provided that [6] one takes the ansatz that Q l assumes its vacuum expectation value 
Ql = y^ dia s(°' 1 '---'V / iV^2,v / iV^T) , (5.13) 

the non-trivial entries of Q 2 are given by (N — 1) complex scalar fields ip and <pi 
(£ — 1, • • • , N — 2) according to 



fe-i = ^' Q^-M- ^ fe €=l,...,JV-2, (5.14) 

and the spatial components of the gauge fields Aj and -B,- (j — 1, 2) are expressed in 
terms of (TV — 1) real-valued vector potentials a = (aj)) and = (bj) (j = 1,2; £ = 
1, • • • , N — 2) satisfying 

A j = B j = diag(0,b?- 2 ,-",b 1 j ,a j ), j = l,2. (5.15) 

Within the above described formalism, the non-Abelian BPS vortex equations 
obtained by Auzzi and Kumar in [6], without restricting to the radially symmetric 
configurations, are of the form 

(d 1 + id 2 )i; = Ka-b 1 )^, (5.16) 
(dx + ityfa = i(b e -b e+1 )fa, l<£<N-3, (5.17) 
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(d 1 +id 2 )<j> N _ 2 = \b N - 2 4> N ^ (5.18) 
a 12 = 2(iV-lV 2 (l-|^| 2 ), (5.19) 
b\ 2 = 2(iV-2y 2 (l + H 2 -2|0 1 | 2 ), (5.20) 
b[ 2 = 2(iV-l-£)a 2 (l + £|0 £ _ 1 | 2 -(£+l)|0 £ | 2 ) ) 

2<£<N-2, (5.21) 

where 

a = (a 1 ,a 2 ) = a 1 + ia 2 , b e = (b[, b l 2 ) = b[ + ib e 2 , £ — 1, - • • ,N — 2, 

are the conveniently complexified gauge vector fields and the summation convention 
is not applied to the repeated index £. As before, the structure of the equations 
(5.16)-(5.18) implies that the zeros of the fields ip,(fie (£ — 1, ■ ■ ■ , N — 2) are discrete 
and of integer multiplicities which may collectively be expressed in the form of the 
respective finite sets 

Zl> = {pi,i,---,Pi,n 1 }, = {pt+i,i,---,Pt+i, nt+1 }, £=1,---,N -2. (5.22) 

For the prescribed sets of zeros given in (5.22), we are to construct a solution of 
(5.16)-(5.21) to realize these zeros. For this problem, here is our main theorem. 

Theorem 5.1. Consider the general BPS system of multiple vortex equations (5.16)- 
(5.21) for (if;,<f>£,a,br) with the prescribed sets of zeros given in (5.22) so that if),(f>£ 
have ni,n e+1 , £ = 1, ■ ■ ■ , N — 2, arbitrarily distributed zeros, respectively. 

(i) For this problem over a doubly periodic domain Q, a solution exists if and only 
if the following (N — 1) conditions 

N-l 



AnJ2 n k < (5.23) 

j N-l N-l \ 

— E^ + aT^E^ < 2X M, (5.24) 

k=l k=2 J 



k=l 

An 



N 

N-l „ N-l _ N-l 



An 



j^Y, nk + jfz^Y, nk + --- + jfr}Y, nk ) 

k=l k=2 k=t J 

£ = 3,-..,JV- 1, (5.25) 



are fulfilled simultaneously. Moreover, whenever a solution exists, it is unique, 
(ii) For this problem over the full plane M. 2 subject to the boundary condition 

H 2 , |0 £ | 2 ^1 as \x\->oo, £=!,■■■ ,N -2, (5.26) 

there exists a unique solution up to gauge transformations so that the boundary be- 
havior stated above is realized exponentially rapidly. 
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In either case, the excited total vortex fluxes are quantized quantities given explic- 
itly by the formulas 



$a= a l2 dx = 2vr^n fe , (5.27) 
J k=i 

/N-l 
b\ 2 dx = 27rJ2 n k, (5-28) 

k=2 

/N-l 
b e 12 dx = 2vr ^ n fc , e = 2,---,N-2, (5.29) 
i — i? 1 1 



respectively. 



To approach the problem, we proceed as before. Firstly, note that, away from the 
zero sets given in (5.22), we may resolve the equations (5.16)— (5.18) to find 

a = -2id(\nip + ln0i + ln0 2 H hln0Ar_ 2 ), (5.30) 

b 1 = -20(ln0 1 + ln0 2 + --- + ln0 JV _ 2 ), (5.31) 

b l = -2i9(ln0 £ + --- + ln0jv_ 2 ), i = 2, • • • , N - 3, (5.32) 

b N-2 = _ 2 id\n<f) N _ 2 . (5.33) 

Thus, following the same procedure as before to substitute (5.30)-(5.33) into (5.19)- 
(5.21), we arrive at 

A(ln|Vf + ln|0 1 | 2 + --- + ln|0 JV _ 2 | 2 ) = X(N - l)(|Vf - 1), (5.34) 
A(ln|^ 1 | 2 + --- + ln|^ 7V _ 2 | 2 ) = A(iV-2)(2|0 1 | 2 -|^| 2 -l), (5.35) 
A(ln|^| 2 + --. + ln|^_ 2 | 2 ) = X(N -1 -£)([£+ 1]|0,| 2 - £|0,_!| 2 - 1), 

£ = 2,---,N -3, (5.36) 
Aln|0 w _ 2 | 2 = A([7V-1]|^_ 2 | 2 -[7V-2]|^_ 3 | 2 -1), 

(5.37) 

away from the zero sets (5.22), where A = 4/z 2 . Next, set m = N — 1 and 

•U! = In |V| 2 , u t = ln|0 £ „x| 2 , e = 2,~-,N-l = m. (5.38) 
Then, the equations (5.34)-(5.37) are converted into 

rn n k 

A( Ml + M2 + --- + Mm ) = Am(e Ul -l)+4vr^^5 Pfcs (x), (5.39) 

k=l s=l 

A(u 2 + --- + u m ) = A(m - 1) (2e U2 - e" 1 - 1) 

+4^J]V lS W, (5-40) 

k=2 s=l 



20 



A(u e + ■■■ + u m ) = \{m-£ + 1) {£e Ue - [£ - lje^" 1 - 1) 

+4^^V iS W, € = 3,...,m-l, (5.41) 

fc=£ s=l 

Aw m = A (me"" 1 - [to, - lje"™- 1 - 1) 

(5.42) 

s=l 

The above formalism may be viewed as an SQCD extension of those of Hong, 
Kim, and Pac [28], Jackiw and Weinberg [30], and Dunne [10, 11], of the Abelian and 
non-Abelian Chern-Simons-Higgs theory. See also the survey [29]. 

In the following two sections, we first consider the equations over a doubly periodic 
domain. Then we consider the equations over the full plane. 

6 Proof of existence for doubly periodic case 

Now consider the equations (5.39)-(5.42) defined over a doubly periodic domain, say 
Vt. We use the direct method to solve them. Let w° be some doubly periodic source 
functions over Q satisfying 



Au2 = -^ + 47r£* w ,»» e=l,2 r ..,m. (6.1) 

I I s=l 

Then the substitutions ue — u® + ve, £ — 1, 2, • • • , to, give us the regularized equations 

. m 

A( Vl +v 2 + --- + v m ) = Am(e u ? +1 ' 1 -l)+-^^n A; , (6.2) 

' ' k=i 

A(v 2 + --- + v m ) = A(m-l)(2e u ^ 2 -e u ? + ^-l)+^f> fc , (6.3) 

' ' k=2 

A m 

A(v e + --- + v m ) = \(m-£ + l) (£e u ° +ve -[£- lje^-i^- 1 - l) + n ^ 

' ' k=i 

j = 3,---,m-l, (6.4) 
Av m = A (me u °™ +Vm - [m - lje"™-^™- 1 - lj + (6.5) 

Integrating the above equations, we obtain the constraints 

[e u " +Vl dx = M-^-l -VnJ = 77x > 0, (6.6) 
Jn \ m k=i J 

C 4 / 1 m 1 m \ 
e u ° +v >dx = |n|--£ -VV + :UU 2 >0, (6.7) 
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/ 

Jn 



e v!i+vt dx 



\n\ - 



4tt / 1 



m ^ m 1 m \ 

— V n fe + + + — — V n fc 

m /L ^' m — 1 m — £ + 1 / 

fc=i fc=2 k=e / 



= r]e>0, £ = 3,---,m, 



(6.8) 



which are the conditions stated in (5.23)-(5.25). Moreover, it will be convenient to 
introduce the transformation 



Wi 

w 2 



v\+v 2 -\ V v m , 

V2 H \~v m , 



v e H N m , 



Vi 
V 2 



Wi - w 2 , 
w 2 - w 3 , 



w e - w e+u 



(6.9) 



ill ^m, i )7) . 

Consequently, the governing equations become 

. m 

Aim = mA (>? + "' 1 - ,i ' 2 - l) + ^n fe , 

' ' fc=i 

Aw 2 = (m - 1) A (2e u 2 +w ™ - e^ 1 - 5 " 2 - l) + ^ n k , 
= (m-£+ 1)A ^e n ° + ^~^+ 1 - [£ - \] e u U+^-i-^ _ i j 

^1 77" 

' ' k=e 

Aw m = A (me< +w ™ - [m - l] e <-^ w ^~ w ^ - l) + y^yn m , 



(6.10) 
(6.11) 



(6.12) 
(6.13) 



whose variational functional is seen to be 



1 



/ n [2mA 
+J(wi,---,w m ), 



Vmil ' + - + 2( m - 1 f+l)A |Vtt " l ' + - + ^ |Vt ^ 1 ' 



dx 
(6.14) 



where 

J(W!,- ■ ■ ,w m ) 



A III' 

1 TTTTT 2^ U k 



1 1 fc=i 



+ ••• + 



47T 



(m-£ + 1)A| 



we 



H — • + me 



47T 

1 - XM nr ' 



w m I > dx. 



(6.15) 
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On the other hand, in view of (6.9), we obtain after some algebra the representa- 
tion ^ 

J(w u ■ ■ ■ ,w m ) = J^l ( jf e""+^ dx - m v)j . (6. 16) 

Thus, we may use the same direct minimization method as before in a verbatim way 
to establish the existence and uniqueness of a critical point of the functional (6.14). 

For completeness, we now sketch how to establish the existence of a critical point 
of (6.14) by a constrained minimization approach. For this purpose, we rewrite (6.6)- 
(6.8) as 

J e (w 1 ,---,w m )= [ e u " +ve dx = r li , £ = !,-■■, m, (6.17) 
Jq 

and consider the problem 

min {l(wi, ■ ■ ■ , w m ) \ (w±, ■ ■ ■ , w m ) satisfies (6.17) and lies in W^ 1,2 (f2)} . (6.18) 

We use the notation w = (wi, ■ ■ • ,w m ). If w is a solution to (6.18), then there 
are some numbers (the Lagrange multipliers) Ai, • • • , X m e M. so that 

(DI(w) + \ 1 DJ 1 (w) + --- + \ m DJ m (w))(w)=0, Vw. (6.19) 

Now insert in (6.19) the test configurations w = = (Su, • • • , S m e), £ = 1, • • • , m. 
We have, after applying (6.17), the relations 

Ai?7i = 0; -A^_i^_i + \ii] e = 0, £ = 2,---,m, (6.20) 

which lead us to Ai = • • • = \ m = 0. In other words, the constraints do not give rise 
to the undesired Lagrange multiplier problem so that a solution of the constrained 
minimization problem (6.18) is a critical point of the functional (6.14) itself. 
Moreover, from the constraints (6.17), we have 

Hi = mi-m+i = hi%-ln ^ e^+^-^dxj , £ = 1, • • • , m - 1, (6.21) 
v m = w m = \ni lm -\n( ( dx \ _ ( 6>22 ) 



n 



Inserting these into (6.16) and applying the condition r)£ > (£ — 1, • • • , m) and 
the Jensen inequality, we again arrive at the coerciveness for the functional (6.14), 



m » 

I(wi, ■ ■ ■ , w m ) > Ci ^2 / \Vw e \ 2 dx-C 2 
i=i J n 



(6.23) 



where Ci,C 2 > are some irrelevant constants. Consequently, the existence of a 
solution to the problem (6.18) follows as before. 



23 



7 Proof of existence for planar case 

To proceed, we define u° e and he as in (4.3) and (4.4) so that £ runs from 1 through m. 
Thus, in view of the translations ui = + v e {£ — 1, • • • , m) and the transformation 
(6.9), the governing equations (5.39)-(5.42) become 

m 

Aw ± = m\ (e u ° +w ^ - 1) + J2 h k, (7-1) 

k=l 

m 

Aw e = (m-£+ 1)A (^ e u ° +w ^ w ^ -[£- i] e «?-i+«*-i-«* - l) + ^ h k, 

k=i 

£ = 2,---,m-l, (7.2) 
Aw m = A (me< +Wm - [m - l] e <-i +w ™- l ~ Wm - l) + h m , (7.3) 

which are the Euler-Lagrange equations of the functional 

I(WU ■■■,Wm) = [ 2 |^E (m _ 1 f+1) l V ^| 2 } ^ + J <^> ->"<»)> 

where 

J(w x , ■ ■ -,w m ) = 

' m—l 

/ „ \ / 

me ra — iu m > ax 



\ \ Yl (£e u " +We ' We+1 - £e u " - w^j + (me" 2 - 

+ / lE 7 T7TT f hk ) ™< + f dx - ( 7 - 5 ) 

i R2 [A^(m-! + l) ^ y A J 



It is important to note that, using the relation 



m—l 



^2 We = ^2 £{wi - w e+ i) + mw m , (7.6) 



e=i 



we can rewrite J(w\, • • • , w m ) defined in (7.5) as 
J(w u ■ ■ -,w m ) = 

J I £ (f° e+Wl ~ Wl+1 - e u? - [w t - w e +i]j + m (e< +Wm - e< - u> m ) J dx 
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Consequently, after some algebraic manipulation, we obtain 



(DI(wi, • • • , w m ))(wi, • • • , w m ) = - 



+ 



(ie u ° [ e "+* - 1] [w e - w e+1 ] + £ [e u " - l] [w e - w e+1 ]j 
+me u ° m [e Wm - 1] w m + m |e u ™ - 

j m— 1 j / m 

+ A^(m-^+l) (E 

£=1 v ; \fc=f 



TO — 1 



M ™* + jh m w m } dx. 



{11 



On the other hand, in view of the transformation (6.9), we have 

m m m 

Ci ^ v] < ^ w £ < c ?Yl v t> 



1=1 

m 



£=1 



ci^|V^| 2 < ^|V^| 2 <c 2 ^|V^| 2 , 



(7.9) 
(7.10) 



i=i i=i i=i 

where ci,C2 > are some constants. Thus (7.8) and (7.10) enable us to arrive at 

(DI(w 1 ,---,w jn ))(w 1 ,---,w m ) > 

m „ m „ 

co / I V ^! 2 + E / { £ ( eU ° e - !] ^ + [ e "° - i] ^) + 9fl>e} dx, (7.11) 

where Co > is a suitable constant and the functions g/s are some linear combinations 
of the functions h^s. It has been seen that the structure of the right-hand side of 
(7.11) indicates that there are constants Ci,C 2 > such that 



(DI(w 1 ,---,w m ))(w 1 ,---,w m )>C i y / / |V^| 2 + -^ dx-C 2 . (7.12 

i { .) -■' \ 1 + \ve\/ 



Therefore, applying (7.9), (7.10), and (7.12), we can again conclude with the coercive 
lower bound 



[DI(wi, ■ ■ ■ ,W m ))(wi, ■ ■ -,W m ) > C 3 y^ ||w<||wi,a(R3) - C 4 , 



(7.13) 



for some constants C3, C4 > 0. Hence, the existence of a critical point of the functional 
(7.4) in the space W^ 1,2 (IR 2 ) follows. Since (7.4) is strictly convex in (w\, ■ ■ ■ ,w m ) G 
VI/ -1 ' 2 (M 2 ) and C 1 , it may have at most one critical point in W^ 1,2 (IR 2 ). 

The rest of the analysis regarding asymptotic estimates and computation of fluxes 
is similar to that of Section 4 and is thus omitted. 
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8 Further applications of direct methods 



In this section, we show that our direct minimization methods may be used to study 
other non-Abelian BPS vortex equations of similar structures arising in SQCD. We 
will present two examples as further illustrations. 



8.1 Vortices in an SO(2N) theory 



In this subsection, we use the direct method developed earlier to strengthen the 
existence results obtained in [35] for an SO(2N) BPS vortex problem formulated in 



, the Lagrangian density 



Recall that, in the work of Gudnason-Jiang-Konishi 
of the non-Abelian Yang-Mills-Higgs model reads 



e 

T 



q}t°q f 



g_ 

2 



2 2 

q}t a q f 



for which the gauge group G is of the general form G = G' x U(l) where G' is 
a compact simple Lie group which may typically be chosen to be G' = SO(2N) 
or G' = USp(2N) (the unitary symplectic group). Assume that a — 1, • • • , dim ((17') 
labels the generators of G', the index indicates the U (1) gauge field, / = 1, • • • , iV flavor 
labels the matter flavors or 'scalar quark' fields, qj, all are assumed to lie in the 
fundamental representation of G'. The gauge fields, gauge-covariant derivatives, and 
field tensors are given by 

A^Al + A^e, V„q f = d M + iAtff, = d„A v - + i[A M , A v ], (8.2) 

respectively, where the generators of G' and U(l), i.e., {t a } and t°, are normalized to 
satisfy 

Tr(t a t b ) = l5 ab , t° = ^^l27v, (8.3) 



/4N 

with l m denoting the m x m identity matrix. When the number of matter flavors is 
^flavor = 2N, the scalar quark fields may be represented as a color-flavor mixed matrix 
q of size 2N x 2N. Restricting to static field configurations which are uniform with 
respect to the spatial coordinate x 3 , a Bogomol'nyi completion [8] may be performed 
to yield the BPS [8, 44] vortex equations [5, 12, 22, 38, 47] 



F? 



12 



'AN 



Viq + iV 2 q 



0. 



(Tr(qql)-vZ) = 0, 



F& a - 9 -{qJ - J\qq ] ) T J) = 0, 



(8.4) 
(8.5) 

(8.6) 
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where J is the standard symplectic matrix 

J =(-°i„ V)- < 87 » 

In its general form, the system of the non-Abelian BPS vortex equations (8.4)- 
(8.6) appears hard to approach and an ansatz-based reduction may be made as a tool 
for further simplification. In the case when G' = SO(2N), the ansatze presented in 
[22] gives us the following matrix forms for the Higgs field, 

/ $l 2JV - 2 \ 
q=\ 0, (8.8) 

\ o ip / 

where <£>, 0, ip are three complex scalar fields, and for the gauge potential, 

Aj = ajl 2N + 6 i diag{0 2 v-2, 1, -1}, j = 1, 2, (8.9) 

where aj and (j = 1,2) are real- valued vector fields. Then, in terms of the com- 
plexified field a and b defined in (2.16), 

the BPS system of vortex equations found in [22] assumes the form 



(9$ = ia$ 


1 


(8.10) 


d<p = i(a - 


f b)4>, 


(8.11) 


dip = i(a • 


-b)ip, 


(8.12) 


e 2 

ai2 = w 


(e 2 -2(iV-l)|$| 2 -|0| 2 -^| 2 ), 


(8.13) 


q 2 

bu = 9 j( 


h/f-H 2 )- 


(8.14) 


Thus, we can recast the system 


of equations (8.10)-(8.14) into 




Aln|$| 4 


= A(ln|0| 2 + ln|Vf), 


(8.15) 


A(ln|0| 2 + ln|^| 2 ) 


= ^(2(iV-l)|$| 2 + |0| 2 + H 2 -a ) 


(8.16) 


A(ln|0| 2 -ln|^| 2 ) 


= 2/(|0| 2 -H 2 ), 


(8.17) 



where we have stayed away from the possible zeros of the fields $, 0, ip. We extend 
our study in [35] and consider a solution so that the zeros of $ coincide with those 
of and ip. As a consequence of the boundary condition, we see that (8.15) leads us 
to the simple relation 

I$| 4 = |0| 2 H 2 . (8.18) 

We are interested in constructing solutions over a doubly periodic domain, Q. The 
multiple vortices are generated from the sets of zeros of and ip, prescribed as 

Z<t> = {pi,---,Pm}, Z f = {gi,- ••,?„}. (8.19) 
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Therefore the vortex-governing equations are then given in terms of the new functions 
u — In |0| 2 and v = In \ip\ 2 as 

m n 

A(u + v) = a (2(N - l)e^ u+v ^ + e u + e v - + AnJ2^)+^J2 6 i^' ( 8 ' 2 °) 

s=l s=l 

n n 

A(u-v) = (3(e u -e v )+4nJ2S P M-^J2 6 ^> ( 8 ' 21 ) 

s=l s=l 

where a,/3,7 are positive constants given by 

« = ^, P = 2g 2 , 7 = e 2 . (8.22) 

In [35], we proved an existence and uniqueness theorem for the solution of (8.20) and 
(8.21) when = (or n = 0) under the necessary and sufficient condition 

Anm (! + !)< 7 |Q|, (8.23) 



and we found suitable conditions under which the solution may be constructed by 
a constrained minimization method. The general existence proof in [35], however, 
is based on a priori estimates and a degree theory argument which is unfortunately 
non-constructive. Here we show that the solution can actually be obtained by the 
(constructive) direct minimization method used in the earlier sections of the present 
paper. For the broadest generality, we consider the presence of the zeros of ip as well 
(n > 0). We are able to obtain the following sharp results. 

Theorem 8.1. For the non-Abelian vortex equations (8.20) and (8.21) defined over 
the doubly periodic domain Q, a solution exists if and only if the condition 

4?r '- + - — - — - <7|fi|. (8.24) 



a (3 

Furthermore, if a solution exists, it is unique and may be constructed by a direct 
minimization method. 

To proceed, let u and v be solutions of the equations 

8=1 ' ' S=l ' ' 

Then the substitutions u = u + U and v — v + V change the equations (8.20) and 
(8.21) into 

A(U + V) = a (2(N - l)e^ Uo+Vo)+1 ^ u+v) + e U0+u + e V0+v - ^ 
An 

+r^(m + n), (8.26) 
Ait 

A{U -V) = (3 (e uo+u - e vo+v ) + — (m-n). (8.27) 
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Next, use the transformation 

f / =\{U + V), l U =f + g, 

{ g = w-n l v =f- 9 . (8 - 28) 

We see that /, g satisfy 
Af = a f(N - l)e^ UQ+VQ)+f + i e Uo+/+s + h Vo+f - g ^r(m + n), (8.29) 

V Z Z Z I i L\ 

Ag = ^( e U0+f+9 -e v « +f - 9 ) +-^-(m-n), (8.30) 
which are the Euler-Lagrange equations of the functional 

I{L9) = J a {^ |V/|2 + 2^ |V ^' 2 + {N ~ 1 ) e4( ' ,o+B0)+/ 

+i (e"°+' + * + e^-«) + + n] - |) / + J^T ["» - n]g\ dx. (8.31) 

On the other hand, integrating the equations (8.29) and (8.30), we obtain the 
natural constraints 

/ \(N -l)e^ uo+vo ^ +f + e U0+f+9 \ dx = l\n\-—( m + n)-^(m-n) 
J n y. ) 2 a [3 

= 77! > 0, (8.32) 

I i(N- l) e ^ u » +v ^ +f + e v » +f - 9 ) dx = ^\n\-—(m + n) + ^-(m-n) 
J n l ) 2 a [3 

= r] 2 >0. (8.33) 
In terms of the quantities 771 and r] 2 in (8.32) and (8.33), we may rewrite (8.31) as 

/(/,«?) = ^{^|V/| 2 + ^|V^| 2 + (iV-l)e^' + ^}da; 



+1 Q { e vo+u+u + e ,oH-v+v| d;r _ ^ _ m v^j . ( 8 . 



34) 



In view of the analysis presented earlier (e.g., §3.1 and (3.17)-(3.22) in particular), 
we see that the existence of a critical point of the functional (8.31) follows as a 
consequence of the condition 771,772 > 0. 
Finally, adding 771 and 772, we have 

V = lM (m + n)>0. (8.35) 

a 

Thus, in view of 771,772 > again, we have 

An 

-77 < j(m-n) < 77, (8.36) 

which leads to the single condition (8.24). On the other hand, it is obvious that (8.24) 
also implies 771,772 > 0. Thus the proof of the theorem follows. 
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8.2 Vortices in a softly broken SQCD model 

In this subsection, we construct multiple vortices in the SQCD model of Marshakov 
and Yung [38] in which the confinement is achieved through Abelian fluxes generated 
in the Cartan subalgebra sector of SU(3) so that, in terms of the Gell-Mann matrices 
A 3 and A 8 , the gauge field assumes the form 

A, = A^X 3 + A^X 8 , (8.37) 

where A^ and A^ are two real- valued vector fields. As in [38], we use u and d to 
denote the up and down colors of quarks, which are represented by a pair of complex- 
valued Higgs scalar fields, say <j)^ and (f>^\ respectively. The u- and d- fluxes will be 
the Cartan subalgebra valued which are induced from the gauge fields 

A W = yl A m + i A W A^ = -^Af + -Af. (8.38) 

With the above notation, the effective action functional of the SQCD model of 
Marshakov and Yung [38] reads 



Q = [ / J_ p(3) F (3)»v , J_ p(8) F (8)iw 

J \ v ^ V 



+ 9 j (i0 (u) i 2 - i0 (d) r - *[i - ^]) 2 + f 4 (i0 (u) r + i0 (d) i 2 - *[i + -]) 2 } dx, 



.39) 



where g,£,u > are all physical constants, and 



V^ d) =^--^ U ' d) , ^ = 0,1,2,3, (8.40) 

are gauge- covariant derivatives to be operated upon (f)\ u ^) , accordingly. Concentrat- 
ing on static case for which the field configurations are uniform in a spatial direction, 
say x 3 , we see that the method of Bogomol'nyi [8] may be used to show that the 



v (u,d) ( u ,d) +iV (u,d) ( U! d) = 

^ ( 2 3) + ^(i0 (u) r-i0 (d) i 2 -ai-^) 

^ 8) + 7/^(l^ (U) | 2 + l^ (d) | 2 -^+-]) 



f 

2V3 

subject to the boundary condition 

(u) | 2 ^e, I0 (d) | 2 ^, 



[8, 44] system [c 


18] 


= 0, 


(8.41) 


= 0, 


(8.42) 


= 0, 


(8.43) 




(8.44) 
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so that the (minimum) vortex-line energy or tension may be calculated via the flux 
formula 

T = 4= [ (F^+uF^)dx. (8.45) 



V3 

With the notation 

= (U) , ^ = (d) , a, = \Af + ^Af, b, = Af\ (8.46) 
and the relations 

4 U) = ^ - \b„ vi u) =d ll - 1 -[a ll - -±=b^ , Vi d) = d, - ±=b„ (8.47) 
we may rewrite the equations (8.41)-(8.43) as 

(d 1 + id 2 ) ( p = i(ja 1 + ia 2 ]--^[& 1 + i& 2 ]^0, (8.48) 

(d 1 + id 2 )^j = ^=(b 1 + ib 2 )^, (8.49) 

a» = ^(e-|0| 2 ), (8.50) 



1 



+ I|0| 2 -|Vf)- (8.51) 



Assume the sets of zeros of and ip are as prescribed in (8.19). Then, as before, 
the substitution, u = In |0| 2 and v = In | | 2 , allows us to transform the equations 
(8.48)-(8.51) into 



i 2 m n 

A M + -At; = |(e«-0+47r^^( a ;) + 27r^^(x), (8.52) 



s=l s=l 

2 



At; = ^(- e u + 2e"-£[2 W -l]) + 47rJ]<5 ? >)- (8.53) 

6 8=1 

We consider a doubly periodic domain Q first. Let u and t> be given in (8.25). 
Then u = u + U and v — v + V recast (8.52)-(8.53) into 

AU+'-AV = ^(e— - + l= + ^, (8.54) 

Ay = £ (-e*° +u + 2e w+v - ?[2w - 1]) + (8.55) 
Set W = C/ + \V. We arrive at 

.2 -1;t// 



A y = |-^' +y -e^ y+H/ -e[2^-l]) + ^, (8.56) 
AW = ^ e u °-5 y+w -e) +^r(2m + n), (8.57) 
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which are the Euler-Lagrange equations of the functional 

I(V,W) = [ |^|VVf + ^|W| 2 + e M °-5 v ' +H/ + e 1 ' 0+v ' 



V 1 g 



i 

U ~2 



g 2 \n\j V # 2 l fi 



(8.58) 

To proceed further, we integrate (8.56) and (8.57) to obtain the constraints 



/ 

Jq 



e u °-^ v+w dx = - %(2m + n) 



9 2 



J 

Jq 



= r?! > 0, (8.59) 
e Vo+v dx = £uM - -J(m + 2n) 

g 2 

r] 2 > 0. (8.60) 



I(V, W)-\ I { ^\VV\ 2 + \WW\ 2 \ dx 



Therefore, we have 

^ - 

9 2 Jn U 

= (^j e U0+LL+iJ dx - Vl lLj + (^f e vo+ ^ +v dx - V2 v}j . (8.61) 

In view of (8.59)-(8.61), we see that the existence of a critical point, in fact, a global 
minimizer, of the functional (8.58) in W 1,2 (Q) follows as before. 

We now consider (8.52) and (8.53) over the full plane. In (8.19), use the notation 

p s =Pi,s, s = 1, • • • ,m = ni; q s = p 2 , s , s = 1, • • • , n = n 2 . (8.62) 

Let v!j and he be defined as in (4.3) and (4.4), respectively, I = 1,2. Introduce the 
translations 

u = ln£ + u1 + vi, v = \n(£w) + u° 2 + v 2 , (8.63) 
and the refined parameters 

« = \g% P = \g 2 i, 7 = |»V (8-64) 

We do so since the ranges of these parameters will not be important for our existence 
theory over M 2 . The equations (8.52) and (8.53) now become 

A Vl + ^Av 2 = a (e"? + ^ - l) + h 1 + h 2 , (8.65) 
Av 2 = -(3 (e u ° +vi - l) + 7 (e u * +V2 - lj + h 2 . (8.66) 
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Set v i + \v2 = wi, v 2 = w 2 . Then we have the following modified system of equations 

A Wl = a (e M ? +,Ul -5 W2 - l) + ^ + -h 2 , (8.67) 
Aw 2 = -p (e u ° +wl ~^ W2 - 1) + 7 (e M "+^ - l) + &2, (8.68) 
which are the Euler-Lagrange equations of the functional 

I( Wl ,w 2 ) = [ ( — |Vwi| 2 + ^-\Vw 2 \ 2 + e u " +Wl ~^ W2 -e< - Wl - -w 2 
J R 2 [2a A/3 I 2 

(8.69) 

It is clear that this functional is C 1 and strictly convex over W 1,2 (WL 2 ). Besides, we 
have 

(DI{w 1 ,w 2 )){w 1 ,w 2 ) = I \-\V Wl \ 2 + ^-\Vw 2 \ 2 \ dx 

+ £ I ( e »S + «- J« _ ^ _ ^ + 7. ( e «i + « _ t ) W2 

+ o~ + 2^ 2 ) Wl + 2^^} dX ' ^ 8 ' 70 ^ 
Thus, we can show that there are constants Ci, C 2 > such that 

(D/(Wl,W 2 ))(Wi,W2) > Cl (||Wl|| W l,2( M 2) + ||w 2 ||wi.2(R2)) - C 2 , WW U W 2 . (8.71) 

Consequently, the existence and uniqueness of a critical point of the functional (8.69) 
over W l,2 (R 2 ) is established. 
In summary, we may state 

Theorem 8.2. Consider the BPS system of SQCD vortex equations (8.48)-(8.51) 
for ((p, ip, aj, bj) with the prescribed sets of zeros given in (8.19) so that 0, ip have m, n 
arbitrarily distributed zeros, respectively. 

(i) For this problem over a doubly periodic domain Q, a solution exists if and only 
if the condition 

max |^(m + 2n), 2m + raj < - ^ (8.72) 

holds. Moreover, whenever a solution exists, it is unique. 

(ii) For this problem over the full plane M, 2 subject to the finite-energy boundary 
condition 

|0| 2 ^£, M 2 -»fw, as \x\-> 00, (8.73) 

there exists a unique solution up to gauge transformations so that the boundary be- 
havior stated above is realized exponentially rapidly. 
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In either case, the excited total vortex fluxes are quantized quantities given explic- 
itly by the formulas 

$ a = J ai2 da; = 2?r(2m + n), $ fe = J b 12 dx = 2V3im, (8.74) 

respectively, and the solutions may be obtained by methods of direct minimization. 

We note that the condition (8.72) is simply a suppressed restatement of the two 
simultaneous conditions (8.59) and (8.60). 

We also note that, applying the relations between the gauge fields a^,6 M and 
we easily obtain the Abelian (or the Cartan subalgebra valued) fluxes 

$ (3) = J F$ dx = 2vr(m - n), $ (8) = J F$ dx = 2v / 3~7r(m + n). (8.75) 

Moreover, the u- and d- fluxes may be expressed by the formulas 

$ (u) = j F (u) = 2^3 7rm) $(d) = J F (d) = 2<s /3 7rn> ( 8J6 ) 

which depend on the winding numbers of the u- and d-Higgs fields (p^ and (j>^\ 
respectively, and give rise to the tension or energy of the vortex-lines 

T = 27r(m + un)(, (8.77) 

according to (8.45). 

We also note that our direct method works easily for the classical BPS Abelian 
Higgs vortex equations [31] defined over doubly periodic domains [60] or formulated 
over compact Riemann surfaces [41, 42]. 
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